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Abstract 

In numerical studies of Gowdy spacetimes evidence has been found 
for the development of localized features ('spikes') involving large gra- 
dients near the singularity. The rigorous mathematical results avail- 
able up to now did not cover this kind of situation. In this work 
we show the existence of large classes of Gowdy spacetimes exhibit- 
ing features of the kind discovered numerically. These spacetimes are 
constructed by applying certain transformations to previously known 
spacetimes without spikes. It is possible to control the behaviour of 
the Kretschmann scalar near the singularity in detail. This curva- 
ture invariant is found to blow up in a way which is non-uniform near 
the spike in some cases. When this happens it demonstrates that the 
spike is a geometrically invariant feature and not an artefact of the 
choice of variables used to parametrize the metric. We also identify 
another class of spikes which are artefacts. The spikes produced by 
our method are compared with the results of numerical and heuristic 
analyses of the same situation. 



1 Introduction 

The question of the nature of spacetime singularities in general relativity is 
far from being answered in general. The encouraging progress which has 
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been made on this problem in the past few years has mainly consisted of 
obtaining detailed information about certain phenomena in more or less spe- 
cial classes of spacetimes. The effects which have been studied in a rigorous 
mathematical way include the decoupling of the evolution at different spa- 
tial points near the singularity (see IQ , fT^j , and references therein) and 



the monotone or oscillatory approach to the singularity, which has been in- 
vestigated in the case of spatially homogeneous spacetimes (p4|, p2|, ||23 



pO| , pTp . In numerical work it has been possible to study both together 
(0) HH) HI)- "^^^ fi^s^ of these effects means that an inhomogeneous model 
can be approximated near the singularity by a family of homogeneous models 
parametrized by the spatial coordinates and this is why the second is of wider 
importance. This choice of topics was stimulated by the pioneering ideas of 
Belinskii, Khalatnikov and Lifshitz (BKL) (see 0] and references therein). 

One class of spacetimes which has turned out to be particularly useful in 
developing ideas about the nature of spacetime singularities is that of Gowdy 
spacetimes. These do not show an infinite number of oscillations as the sin- 
gularity is approached in the way that generic spacetimes are supposed to 
do according to the BKL picture. Instead they appear to show a monotone 



behaviour near the singularity which is called velocity dominated ([0], 0])- 
Thus they provide a relatively simple arena for studying spacetime singular- 
ities. On the other hand, as will become clear in the following, the Gowdy 
spacetimes do display quite complicated dynamical behaviour near their ini- 
tial singularities and it is probable that what we can learn from them will be 
useful in understanding more general classes of spacetimes. 

The Gowdy spacetimes which have been studied most extensively are 
those with spatial topology and only these will be discussed in the fol- 
lowing. Grubisic and Moncrief derived a consistent formal expansion for 
Gowdy spacetimes near their singularities under a certain assumption. For 
Gowdy spacetimes it is possible to consider a quantity v called the velocity. 
(Since f is a non-negative scalar quantity it would be more appropriate to 
call it speed, but we will retain the usual terminology.) If velocity domi- 
nance holds this quantity has a (space-dependent) limit as the singularity 
is approached. It is also possible to associate a sign to each spatial point 
and thus define a real-valued function k with the given magnitude and sign. 
This function will be called the asymptotic velocity. The Grubisic-Moncrief 
expansion is consistent provided < /c < 1, which defines what is called 
the low velocity regime. For high velocity {k > 1) this is no longer true in 
general. (Some remarks on the meaning of A; < will be made later.) It was 
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suggested in |T3 that in those cases where the expansion is not consistent 



the solution evolves in such a way that the velocity is eventually less than 
one close to the singularity. It was also noticed that under an additional 
non-generic condition it is possible for the expansion to be consistent with 
high velocity. 

A little later Berger and Moncrief p observed the development of large 
spatial derivatives near the singularity in numerical calculations of Gowdy 
spacetimes. They called these 'spiky features'. We will abbreviate this to 
'spikes'. They presented evidence that these were genuine features of the 
spacetimes and not numerical artifacts. Spikes were seen to be associated 



with the situation where the non-generic condition of |]T4| is satisfied at an 
isolated point. In this case it can happen that the velocity stays high at 
the special point as the singularity is approached while it is forced below 
one everywhere else. This gives an intuitive picture of how under certain 
circumstances strongly spatially inhomogeneous behaviour may develop in 
Gowdy spacetimes as the singularity is approached. 

In the special case of the polarized Gowdy spacetimes velocity dominance 
has been proved [|l6l . There is also a theorem on the structure of the initial 



singularity in those non-polarized vacuum models evolving from data which 
are close to those for a particular Kasner solution (with exponents |, f, — |) 
| |12| |. The model solution has velocity v identically zero. The discovery of 
spikes in the evolution of general Gowdy models implies limits to the way in 
which these results for particular subclasses could generalize. 

A complementary development to these theorems on the nature of the 
singularity arising from given Cauchy data for a Gowdy spacetime is the 
construction of large classes of solutions with singularities of a particular 
type. This was carried out for the case of velocity dominated singularities 
with asymptotic velocity satisfying < A; < 1 by Kichenassamy and Rendall 



17| . Solutions were constructed depending on four free functions, which is 
the same number as in the general solution. Intuitively this suggests that 
solutions corresponding to an open set of initial data are being obtained by 
this procedure, although it has not been possible to prove this yet. Solutions 
depending on three free functions were also constructed under the weaker 



restriction k > 0. In ||T^ the free functions were required to be analytic. 



but this undesirable restriction was removed in |jT9[ . In the following the low 



velocity solutions of |jT3 and [|T9l will be used as a starting point to produce 
another family of solutions, depending on the same number of free functions, 
which exhibit spikes. 
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If one looks at numerically produced spikes (see e.g. figure 2 in |^) then 
it becomes clear that there are two different types and that within any one of 
these types the individual spikes show a great similarity among themselves. 
This is a strong motivation for finding a simple explanation, and the main 
result of this paper is to do this in a mathematically rigorous way. The prin- 
cipal unknowns in the Einstein equations for Gowdy metrics are real- valued 
functions P and Q. In one case P has a spike and Q has a kind of disconti- 
nuity while in the other case P has a spike of the same shape but pointing in 
the opposite direction while Q is smooth. For reasons to be explained later, 
these will be called 'false' and 'true' spikes, respectively. A key feature is 
that in the first case geometric quantities such as curvature invariants show 
no irregular behaviour at the spike as the singularity is approached, while in 
the second case they do. 

The field equations in Gowdy spacetimes are known to be closely re- 
lated to the geometry of the hyperbolic plane (see, for instance, [0). The 



interpretation in terms of hyperbolic geometry leads rather easily to a trans- 
formation (inversion) which produces solutions with false spikes from the 
smooth low velocity solutions of [1^. More subtle is a further transformation 
(Gowdy-to- Ernst transformation) which converts solutions with false spikes 
into solutions with true spikes. It is related to the Kramer- Neugebauer trans- 
formation ([|TB|, []TU|) occurring in the theory of stationary and axisymmetric 
spacetimes. 

The paper is organized as follows. In section 2 the basic notations and 
equations are introduced and the two transformations of key importance in 
the following are defined. The third section introduces the concepts of false 
and true spikes, shows how solutions displaying features of this kind can be 
produced and gives asymptotic expansions for these solutions near the spikes. 
In the fourth section the analytical and numerical results are compared. Sec- 
tion 5 relates what we have found here to the method of consistent potentials 
and hence to the dynamics which produces these solutions. Finally, section 
6 discusses some further developments of the results. 



2 The Gowdy spacetimes 

The Gowdy spacetimes on are solutions of the vacuum Einstein equations 
which can be written in the form 
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They admit U{1) x U{1) as an isometry group, acting by translation of 
the periodic coordinates y and z. The notation here differs somewhat from 
that used in ||T^. The quantities {P,Q,t) occurring here are related to the 



quantities {X, Z, t) of by 



t = e~^ P = -Z Q = X 

The coordinate t, and hence the coodinate r, is determined in an invariant 
way by the geometry: t is proportional to the area of the the orbits of the 
symmetry group. The functions P and Q depend on r and x. For a Gowdy 
spacetime on they should be periodic in x, but most of this paper is con- 
cerned with constructions which are local in x. Information about the initial 
singularity, which is at f = 0, is encoded here in the asymptotic behaviour 
of the metric as r oo. The vacuum Einstein equations imply that the 
functions P and Q satisfy 

Prr = e'^Ql + e-'^iP^^-e'^Ql) 

Qrr = -2PrQr + e'^^{Q,, + 2P,Q,) (2) 

and that A satisfies 

A. = -iP^ + e'^'QD-e-'^iP'^+e'^Ql) 

A, = -2{PrP, + e^''QM (3) 

These equations determine A up to a constant. When the above equations 
are satisfied P, Q and A define a solution of the vacuum Einstein equations 
via (|1D. 

Suppose that a solution {Pq, Qq) of (||) is given. Define functions (Fi, Qi) = 
IiPo,Qo) by 



e--^ = = ^° (4) 



Then also satisfy (0), as can be checked by direct calculation. The 

transformation (FqiQo) —-^ [PiiQi] has been denoted by /, since it has a 
geometric interpretation in terms of an inversion in the hyperbolic plane. 
This will now be explained in some detail. 

The description of Gowdy spacetimes in terms of the three real-valued 
functions P, Q and A of r and x can be replaced by one with greater invari- 
ance which leads to useful insights. For this purpose it is useful to think of 
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these as functions of t and x which, with a shght abuse of notation, will be 
denoted by the same letters. Let M be the manifold with coordinates t and 
X, which is just an open subset of R^. Let the flat metric —dt"^ + dx^ on M 
be denoted by /. Define a mapping from M to by the pair U = (P, Q). 
Let (A^, h) be R^ with the metric h defined by dP^ + e^^dQ"^. The derivative 
dU is a section of the tensor product of the cotangent bundle of M with 
the pull-back via U of the tangent bundle of to M. There is a natural 
connection defined on this bundle by means of the Levi-Civita connections 
associated to the metrics / and h. This defines a covariant derivative, call it 
D. The Gowdy equations for P and Q take the form ii f{DdU) — t^^Ut = 0. 
This a slight generalization of the equation for a wave map (also known as a 
hyperbolic harmonic map or nonlinear cx-model) which is ti f{DdU) = 0. The 
functions of P and Q occurring in the equations for A are also expressible in 
terms of h. They are just h{Ut, Ut) + h{Ux, Ux) and h{Ut, Ux)- It follows that 
all the Gowdy field equations can be expressed in a way which is covariant 
with respect to {N,h). An advantage of this is that applying any isometry 
of /i to a solution of the Gowdy equations will give a new solution. It can be 
shown that the metric corresponding to the functions (P, Q, A) after trans- 
formation is isometric to the original one, the transformation being effected 
by a linear change of the variables y and z. This will not be proved here, 
since we only need it in one special case (inversion), where it can be checked 
directly. The fact that the two metrics are isometric means that quantities 
such as invariants of the spacetime curvature are identical for the original 
solution and the transformed solution. 

Another way of looking at the invariance properties of the equations and 
their relation to wave maps is via a Lagrangian formulation. A wave map 
from two-dimensional Minkowski space to the hyperbolic plane is the Euler- 
Lagrange equation corresponding to the action 



Pd,U%U''hABdV{f) 

= J{{-P^ + P'x) + e'^{-Q',+Ql))dtdx (5) 

The Gowdy equations are the Euler-Lagrange equations for the analogous 
action 

/ ii-P^ + e'^^Pl) + e'^{-Ql + e-'^Ql))dTdx (6) 
One isometry of the metric h is given by the transformation (^) . This was 



written down in a slightly different notation in [17|. Its geometrical origin is 
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as follows. Define Y = e~^. The metric h is that of hyperbolic space and 
the coordinates {Q, Y) give the well-known upper half-space representation 
of that space on the region y > 0. In that representation the transformation 
introduced above becomes inversion in the unit circle and it interchanges the 
origin with the point at infinity. For this reason this transformation is referred 
to as inversion / in the following. It is induced by the linear transformation 
of coordinates which interchanges y and z. Another representation of the 
hyperbolic space is given by the Poincare disk model, where the points of the 
boundary of the half-space model are all on an equal footing with the point 
at infinity. From this point of view it is clear that the coordinates {P,Q) 
break the symmetry of the hyperbolic plane by singling out a point of its 
ideal boundary as the point at infinity. The inversion / changes this choice. 

Next another transformation which generates new solutions from old will 
be introduced. We know of no simple geometric interpretation such as that 
given for / above. In contrast to /, which transforms one coordinate rep- 
resentation of a metric into another representation of the same metric, the 
transformation which will be considered now produces a genuinely new met- 
ric (i.e. one which is not isometric to the original one). This transformation 
was suggested by the fact, already discussed in detail in [|TT], that a given 



Gowdy spacetime defines two pairs of functions which satisfy identical equa- 
tions. These come from the Gowdy and Ernst methods of reduction with 
respect to the Killing vectors of the spacetime and so we will call this trans- 
formation GE for Gowdy-to- Ernst. If {Pi,Qi) is a solution of the Gowdy 

GE 

equations then the transformation {Pi, Qi) — > {P2, Q2) is defined by 

P2 = -Pi + r Q2r = -e'^^'^-^^Qi. Q2X = -e'^'^QiT (7) 

In fact, Q2 is only defined up to an additive constant by these relations and 
it is necessary to be careful about specifying the constant when applying 
the transformation for a specific purpose. The integrability condition only 
assures the existence of a solution locally in x and in most of this paper 
we will restrict to that situation. See, however, the comments on spikes in 
(global) Gowdy spacetimes in section ^. The transformation GE is such that 
if (Fi, Qi) satisfy the Gowdy equations (0) then (F2, Q2) do so too. That the 
integrability condition for Q2T and Q2X is satisfied can be verified directly 
by calculating the x-derivative of the right hand side of the expression for 
Q2r and the r-derivative of the right hand side of the expression for Q2X and 
seeing that they agree if Pi and Qi satisfy the Gowdy equations. There is a 
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close relationship between this transformation and the Kramer-Neugebauer 



transformation for stationary axisymmetric spacetimes (see |]T8[, [|T0[)- They 
differ by no more than signs related to the different character of the Killing 
vectors (both spacelike in the one case, one timelike and one spacelike in the 
other case). 

Since the relations (|^ only fix Q2 up to a constant, there is in fact a 
whole one-parameter family of solutions produced by this procedure which 
differ by translations in the Q-direction. On the other hand, changing Qi by 
a constant does not change Q2- 



3 Construction of solutions with spikes 

The central procedure of this paper is the manufacture of solutions of the 
Gowdy equations with spikes out of solutions without spikes by means of 
the transformations / and GE introduced in the last section. The source of 
solutions without spikes which we draw on is a class of solutions constructed 
in |jl9[ and so we begin by describing those solutions in terms of the variables 
used here. Let e be a positive constant and let k{x), Xq{x), (t>{x) and i/j^x) 
be periodic C°° functions, i.e. smooth functions defined on S^. Suppose that 
e < 2k{x), e < 2(1 — k{x)) and < k{x) < 1. Then there exist unique C°° 
functions u{t,x) and v{r,x), periodic in x, which converge to zero in the 
limit that r ^ 00 such that the functions P{t,x) and Q{t,x), 

p = kr-^-e-'^u Q = Xo + e-'^''^ {ip + v) (8) 

satisfy the equations (0). A similar statement holds locally. If the input func- 
tions k, Xq, (p and ip are defined for x belonging to some open interval, then a 
unique solution as above is obtained for x belonging to some slightly smaller 
open interval. The function k coincides with the asymptotic velocity men- 
tioned in the introduction and is positive in this case. Note that statements 
of this kind were originally proved in [ITTl under the additional hypothesis of 



analyticity of the input functions. We will refer to these solutions (whether 
analytic or not) as the low velocity KR (Kichenassamy-Rendall) solutions. 
It follows from what has been said about P and Q above that 



P = kr-ip + o{l) Q = Xo + e-""^ ij + oie-'"^) (9) 

It is also important for the following to know that these asymptotic expan- 
sions can be differentiated term by term in an appropriate sense. This is a 
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consequence of the fact that the functions u and v are regular in the sense 
of Definition 1 in fTOll. It can be concluded that 



diP = dikr - di^ + o{l) (10) 

for all j > 0. For Q we have the slightly more complicated statement that 

diQ = diXo + diie-''^ ^) + o{T^e-^'^) (11) 

In this relation it is also possible to replace di{e~'^^'^ifj) by dik{—2Ty e~'^'^'^ ip . 
Estimates for time derivatives can be obtained by substituting these expres- 
sions into equations (||). (To be precise, this is true for time derivatives of 
order at least two. The statements for the first order time derivatives are 



direct consequences of the theorems of |T^, which were proved using a first 
order formulation of the equations, where the time derivatives of the basic 
variables are included as unknowns.) There results the estimate 

dldiP = d'Mkr) + 0(1) (12) 

for all J > and / > 0. Of course the first term on the right hand side 
vanishes as soon as / > 1. Similarly 

d[diQ = d[di{e-^^^ ij) + o{r^e-^'^) (13) 

Consider now a low velocity KR solution {Pq,Qq) such that the corre- 
sponding function Xq vanishes at an isolated point x spike- Applying the 
transformation / produces a new solution (Pi,(5i) of the Gowdy equations 
according to equations @). The intuitive idea is that the new solution has a 
spike at Xspike, a statement which will be made precise later. Applying the 
transformation / to (|^) with {P,Q) replaced by {Po,Qo), one finds that 

The aim now is to investigate the qualitative behaviour of this solution in 
the limit r — oo. This will be done by fixing a point x and expanding in 
r in a suitable way. It turns out that the result is quite different depending 
whether x is equal to x spike or not. At x spike one finds that 

= -kr + (f + e~''''ua 
Qi = e-'''{^ + v) + e-'^Va (15) 
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for some functions Ua and Va with hm^_,oo Ua — and HuIt-^^ Va — 0. Else- 



where one finds that 



Pi = kr - (p + IuXq + e~^''ub 

= -^ + e-"'^{-i,-^ + Vb) (16) 



for some functions iih and Vh with hmT-_^oo "^6 = and hm^^oo Vb = 0. At each 
spatial point the asymptotic behaviour in time is of a similar form to that of 
the KR solutions. However the coefficients of the leading terms behave in a 
very non-uniform way as a function of x close to the spike. The coefficient 
corresponding to k{x), the asymptotic velocity, has a discontinuity at Xspike 
while the coefficient corresponding to Xq is unbounded near that point. This 
means that the spatial derivatives of the functions P and Q close to the point 
Xspike grow much faster as r — > cxd than is the case for the KR solutions. This 
behaviour shows up as spikes in numerical calculations. Since it is obviously 
an effect of the parametrization of the metric and not a geometrical one, we 
refer to this kind of point as a 'false spike'. Note that at a false spike the 
asymptotic velocity is negative. This shows that the sign of k has no absolute 
geometric significance in general. 

Applying the transformation GE to (Pi, Qi) gives a new solution (P2, Q2)- 
It also has a spike at the point Xgpike and, as will be proved later, curvature 
invariants show highly non-uniform behaviour there. This shows in particular 
that in this case there has been a change in the geometry and so we refer to 
this kind of point as a 'true spike'. In this case one finds that 

P2 = {l-k)T + ip + e-'^u-\n{[Xo + e-^''^{ijj + v)]^ + e-^''^e^^e^''''^''} 
Q2r = e-^^{-Xo,,- e-^^^{^|J^ + v^) + 

[e^'^Xo. + ij,, + v,,- 2k^T{^ + v)][Xo + e-2fc-(V; + v)fe-''^e-''-'^'' 
-2t,.rXo + 2[Xo + e-^'^{^ + v)][ip^ + e-'^u^]} 
Q2. = 'e--''^Vr + i-2k{t^ + v) + v^)[Xo + e-^''^itlj + v)fe-^'^e~^'-'^'' 

-2kXo + 2Xoe-'^{ur - eu) + 2e-^^''+'^^ {^j + v){ur - eu) (17) 
At Xgpike one finds that 

P2 = {l + k)T-(p + e-'^Uc 
Q2r = -e-^^Xo, + e-^^'^'^^{-ilj, + Xo,ilj^e-^^ + 2ijip, + v,) 
= -e-^''^Vr + 2e-^^'+'^^{i; + v){ur-eu) 

-2ke-^''^iP^e-^'^ + e-^''^Vc (18) 
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with limT-_>oo Uc — 0, lim^_>oo Vc — and lim^_^oo Vc — 0- Elsewhere one finds 
that 

P2 = (1 - k)T + ip- In X^ + e-'^Ud 

+X^{i^^ - 2hTij)e-^'^ - 2k^TXo + 2Xocp^ + Vd} (19) 
Qsr. = -2kXo + {-2k{ij + v)+Vr)X^e-'"^e-^^~'^'' 

-4ke-^''^xjj^Xoe-^'^ + 2X06''^ {ur - eu) + e'^'^^Vd (20) 

with hniT-^oo Ud — 0, \imr^^ Vd — and \imr^^ Vd — 0. 

While the pointwisc expansions which have been obtained so far contain 
much valuable information, it is also useful to derive uniform expansions of 
the solutions on suitable regions. This will only be done for the true spikes, 
which are the main focus of interest here. For a true spike Q2 behaves in a 
regular way and so the point Xspike does not have to be considered separately. 
It is possible to get a uniform asymptotic expansion for r large and a; in a 
neighbourhood of Xgpike- The expansions for P and Q and their derivatives 
imply corresponding expansions for Q2T, Q2X and their derivatives. An ex- 
pansion for Q2 itself can be obtained straightforwardly by integration. All 
these expansions are of the same form as those of the starting low velocity 
solution. 

In the case of P2 it is necessary to distinguish between two different 
regions in order to get uniform expansions. Let k^ be a positive number. 
The interesting case is where k^ is shghtly greater or slightly less than the 
value of k at x spike, and hence has the same property for x close to x spike- 
Let (t(t, x) = e'^*'^|x — Xspikel- Intuitively the function a is a measure of the 
relative influence of the spike and the part of the boundary with low velocity 
asymptotics. We will see that the spike dominates when a is small while the 
low velocity behaviour dominates when a is large. 

For convenience of notation, let Sx — \x — Xspike\- Consider a region Ri 
where r is sufficiently large, 5x sufficiently small and a > 1. We also make 
the genericity assumption that |Xox(:i'6pifce)| = > 0. Given < a < 1 it can 
be assumed by choosing i?i to make 6x sufficiently small, that on the region 
Ri we have |Xo| > arjSx and 

l^or' < a-'rj-'e''*^ (21) 

Now choose A:* < k{xspike)- Then 

[Xo + e-2'=-(^ + v)]^ + e-^'^-e^^e^^"''" = ^^[l + 0(e-2('=-'=*)^)] (22) 
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The logarithm of this expression is then equal to the logarithm of Xq up to 
a remainder of which is 0{e~'^^''~^*^'^). In this way it can be seen that P2 
has an expansion on the region Ri which is of the same form as that for the 
starting function P. In particular this expansion can be differentiated term 
by term. 

Next consider a region R2 where r is sufficiently large, 6 is sufficiently 
small and a < 1. Then we can carry out a very similar procedure. This time 
we choose an arbitrary a > 1 and fc* > k{xspike)- It is convenient to rewrite 
the expression for P2 as 

P2 = (1 + A;)r - (/? - e-'^u - ln{l + [e'^^Xo + e^'="(^ + v)]h-^^e-^^"^''} (23) 

On the region R2 an estimate of the form e'^^|Xo| < ar/e*^'^"'^*-*'^ holds. Thus 
the logarithm in the last expression for P2 is ©(e^^'^"'^*-*'^) on R2. We obtain 
a uniform asymptotic expansion for the solution on the region R2 which can 
be differentiated term by term. This expansion is of the same form as those 
for high velocity KR solutions. 

The significance of the above observations concerning uniform asymptotic 
expansions is that in order to determine the behaviour of curvature invariants 
in the regions Ri and R2 it is sufficient to use the expressions already derived 
in jr^, since all the assumptions needed to do that computation are also 
satisfied on the regions Ri and R2 in the present case, although for rather 
different reasons. In particular the information about what happens as a 
point of the boundary is approached with x a constant other than Xgpike can 
be read off from the asymptotics in the region i?i while the information about 
what happens when the boundary is approached along x = Xgpike is given by 
the asymptotics in the region R2. It follows that the Kretschmann scalar 
K = R°'f^'^^ Ra/3'ys is given up to terms of lower order by K = '^^^B^k) 
where B{k) = 64(A; — + + 3)~^ and k is the asymptotic velocity 

at the given value of x. For x 7^ Xgpike this is k = 1 — k{x) while k{xspike) = 
1 + k{xspike)- For the solutions under consideration B{k) never vanishes. 
Thus the Kretschmann scalar blows up like a more negative power of the 
Gowdy time t at the spike than elsewhere. Since the Gowdy time is fixed 
invariantly by the spacetime geometry, this provides clear evidence that the 
non- uniformity near a true spike is a feature of the geometry and cannot 
be removed by a reparametrization as can be done in the case of the false 
spike. Note that there is a uniform lower bound for the rate of blow-up of K 
in a neighbourhood of the spike and so there can be no violation of cosmic 
censorship there. 
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The mean curvature blows up like This shows that the singu- 

larity is a crushing singularity even in the presence of a true spike. If the 
Kretschmann scalar is expressed in terms of the mean curvature then the 
same power is obtained everywhere. However the constant of proportionality 
changes discontinuously as a function of x as the spike is passed. 

In [|l^] Hern investigated the behaviour of curvature invariants in dynam- 
ical numerical calculations of Gowdy spacetimes. Our results are consistent 
with his and in some cases help to explain them. He investigated more invari- 
ants and more aspects of their behaviour than we have done. The approach 
used here should also suffice to obtain more detailed information about more 
curvature invariants in the spikes we manufacture, if desired. 

4 Comparison with numerical results 

In the previous section we produced solutions of the Gowdy equations where 
the functions P and Q behave in an inhomogeneous way near the initial sin- 
gularity. The motivation for this was to provide an explanation and rigorous 
confirmation of various features which had been observed in numerical calcu- 
lations. In this section the analytical and numerical results will be compared 
so as to show their similarities. 

As already mentioned in the introduction, looking directly at the results of 
the numerical calculations leads to the conclusion that near the singularity, 
and for generic initial data, there are a finite number of points where P 
has a spike and that these points can be classified into two types according 
to whether the spike in P points upwards or downwards and whether Q is 
smooth or also exhibits inhomogeneous behaviour. We claim that the 'false' 
and 'true' spikes constructed in section 3 reproduce the qualitative features 
observed numerically. 

Consider first false spikes. The coefficient of the leading term in the ex- 
pansion of the function Pi has a discontinuous limit at the singularity. The 
asymptotic velocity jumps. At the spike it is between and —1. At fixed 
time we see the infiuence of the expression 2Xq a little away from the spike. 
Under the generic assumption that Xqx{x spike) 7^ this behaves like a con- 
stant times log I a; — Xspike\ and accounts for the shape of the fianks of the 
spike. The spike points downwards. The coefficient of the leading term in 
the expansion for Qi is unbounded near the spike and at fixed time behaves 
like a constant times {x — Xspike)~^ (assuming the genericity condition). This 
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Figure 1: Pi in a neighbourhood of x spike at small t. 
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Figure 2: Qi in a neighbourhood of Xgpike at small t. 

leads to a growing jump in Q between positive and negative values as the sin- 
gularity is approached. Another quantity which has been studied numerically 
is the momentum conjugate to Q in a Hamiltonian formulation, ttq = e^^Qr- 
Calculation of tiqi from Pi and Qi shows that it tends to zero as t — >■ at 
X spike- These features are observed in numerical calculations of the dynam- 
ics. To compare the two, we present pictures of the results of applying the 
transformation / to a low velocity KR solution in figures For purpose of 
illustration we choose periodic functions (with period 27r) for /c, Xq, ip and -i/', 
and ignore the decaying unknowns u and v. The salient features of the spikes 
do not depend on the particular form of the data, but to obtain a spike Xq 
must vanish at an isolated point. 

Consider next true spikes. Applying the GE transformation to the false 
spike shown in figures we get a true spike. Since P2 is obtained from 
Pi by reflection (about the x-axis) and translation (in the vertical direction) 
it is obvious that the profile has the same shape as that of the false spike 
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Figure 3: ttqi in a neighbourhood of Xgpike at small t. 
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Figure 4: ttq^ in a neighbourhood of x spike at small t. 
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shown in figure |I|, except that in this case the spike points upwards instead 
of downwards. On the other hand Q2 tends to a continuous hmit as t — 
and its spatial derivative tends to zero at Xgpike- Since Q2X = — ttqj, its graph 
is obtained from figure ^ by reflection about the x-axis. The most dramatic 
behaviour at a true spike is exhibited by ttq, as can be seen in figure The 
shape of ttq at a true spike can be understood from the identity = —Qix- 
Figures along with P2 and Q2X as obtained from figures |I| and ^, should 
be compared with figure 2 in and figures 2 and 3 in 0. 



5 Comparison with the method of consistent 
potentials 

Numerical simulations of Gowdy spacetimes consistently show that the 
speed of the solution in hyperbolic space, 



P2 + e^PQl (24) 

is, asymptotically, less than one almost everywhere Furthermore, the 
numerical simulations show that w < 1 is obtained through a series of steps. 
On each step v is essentially constant in time. From one step to the next, v 
goes to approximately |2 — i^l, until f < 1 on the final step. In the absence of 
spikes, on each step the numerical solution is approximated well by a geodesic 
loop solution |M, in which the solution follows a geodesic in hyperbolic space 



at each spatial point. A geodesic loop solution is a one parameter family of 
Kasner solutions. Thus the steps are known as Kasner Epochs. 

In this setting the one parameter family of Kasners occurs as a solution 
to the following differential equations for P and Q. 

Ptt 6 Qr 

Qrr = -2PrQr (25) 

The numerical simulations show that, during a Kasner Epoch, the terms on 
the right hand side of the evolution equations (2) which are absent from the 
right hand side of equations (E^p are small. However, unless f < 1 or some 



non-generic condition is satisfied, one of the absent terms eventually grows. 
On the other hand, if v < 1 the numerical simulations show that the absent 
terms are all decaying. The Kasner solutions can be written as follows, with 
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constants a > 0, P > 0, ( and ^ 0. 

P = lna + ln(e-'^" + C^e^") 

A one parameter family of Kasners is obtained by allowing a, f3, C and C, to 
depend on x. Calculation shows that the speed in hyperbolic space for (P, Q) 
as in equations (^) isv = p. If (^), with v < 1 and C 7^ 0, is substituted into 
the right hand side of equations then the neglected terms are all decaying 
for large enought r. (The result in takes into account corrections to all 



orders and finds them to be decaying as r — > cxd.) On the other hand, if (pq), 
with V > 1 and C 7^ 0, is substituted into the right hand side of equations (H) 
then one of the terms in equations (0) that is absent from equations (pS]) 
eventually grows, in agreement with the numerical simulations. 

As the term grows, the one parameter family of Kasners stops being a 
good approximation. The numerical solution is next approximated well by 
a "curvature bounce solution." The name refers to the solution's role as the 
approximation of the transition from one Kasner Epoch to the next, in which 
a term related to the spatial curvature grows and then decays again. In a 
curvature bounce solution, f ^ |2 — t;| 0, in agreement with the numerical 
observations. For a period of time, both the one parameter family of Kasner 
solutions and the curvature bounce solution approximate the Gowdy solution. 
That is, the curvature bounce solution approximates the latter part of a step, 
the transition and the first part of the next step. 

It turns out that the curvature bounce solutions are related to the Kasner 
solutions. Application of the transformation GE to a solution of equations (2) 

GE 

during a curvature bounce, (P, Q) — > (p, q), and then neglect of small terms 



implies that p and q satisfy equations Therefore the expressions on the 
right hand side of equations ( p6D can be used for p and q. The transformation 
GE is its own inverse, so transforming {p,q) back to {P,Q), and defining 
'^Q = -Qx gives 

P = -lna-ln(e-(^+i)^ + CV^-i)^) 
Qx = 2a/3C (27) 



a{l + Ce^f^^) [a 1 + Ce^^^ 1 + Ce^^^ 
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Since (p, q) do not satisfy the full Gowdy equations, the integrability con- 
dition for Q is not satisfied after the transformation back. If it were, then 
it would be the case that e~^^7rQ = Qr- The conditions which must hold 
so that (p, q) approximately satisfy the full Gowdy equations (0) imply that 
the integrability condition is approximately satisfied. Upon substitution of 
equations (pTj) into equations (|^), one finds here too, that if (3 < 1 and C 7^ 0, 
the neglected terms are all decaying for large enough r, but if /5 > 1 and 
C 7^ 0, one of the neglected terms in the equations (0) eventually grows. In 
this case, the next stage of the evolution can then be approximated by a one 
parameter family of Kasner solutions (pB]) and the cycle continues. 

As mentioned in section |^, if no terms are neglected after the transfor- 
mation (P, Q) — > (p, q), then (p, q) satisfy the Gowdy equations (H) exactly. 
In this case the result in [|T^ can be applied to (p, q). Transformation with 
GE to (P, Q) then shows that corrections to (^Tf) to all orders are decaying 
as T CX3, if /? < 1 and C 7^ 0. 

The occurrence of spikes can be understood within this approximation. 
A "downward pointing" spike first appears during a Kasner Epoch if ^ = 
at an isolated point. An "upward pointing" spike first appears during a 
curvature bounce, again if C = at an isolated point. The resulting shape of 
P and Q as given in equations (|26|), and of P, Q and ttq in equations (pTj), 
qualitatively agree with the numerical observations. If the approximation 
remained a good one, this would imply that the range of A; at a spike is 
greater than obtained in section ^ 

To consider the magnitude of k, let Pn = Pq + ■ ■ ■ + Pn and Qn = 
Qo + ■ ■ ■ + Qn- Let (Fo, Qo) be given by equations (|26| ) with C = at xq- At 
large r, Fqt ~ —P at Xq. Let Pn+i be obtained by substituting Pn+i into 
the left hand side of equations (0) and {Pn, Qn) into the right hand side and 
Qn+i be obtained in the analogous way. 

At Xq, Pirr is decaying for large enough r if < /3 < 1 or for any j3 > 
if the additional condition, = at xq, is imposed. However, if /? > 1 and 

7^ at Xq, then Pi^-j. is growing exponentially for large r and is strictly 
positive. This suggests that k < —1 will not occur in general at a spike. 

Calculation of P2tt shows that it is growing exponentially at Xq for large 
T a P > 2, = at Xq and Cxx 7^ there. The dominating term is, again, 
positive definite. If (xx does vanish at xq and if /? < 3, then P2rT is decaying 
for large enough r. Inspection of the differential equation satisfied by Pnrr 
suggests that this pattern will continue. In the expression for Putt there is a 
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term of the form +(d^CYe^^^~"'^'^ . If cancellation of this term does not occur, 
this suggests that k < —n will not occur at a spike without conditions on 
the spatial derivatives of C up to order n. 

The conditions on the spatial derivatives of C can be expressed as condi- 
tions on ttq. Let ttqo = e^^^Qor- From (^) it follows that ttqo = — 2a/?C- If 
C = 0, then d'^TiQQ = is equivalent to = 0. The analysis suggests that 
k < —n does not occur at a spike unless ding vanishes asymptotically at the 
spike, for all integers < j < n. 

Because of the relationship of equations ( P7| ) to equations (^), the analy- 
sis of a spike which begins when the evolution is approximated by a curvature 
bounce solution can be obtained from the above by replacing the conditions 
on the spatial derivatives of vtq with conditions on the spatial derivatives of 
Qx, and replacing P.^ ~ — at large r with P^- ^ 1 + f]. So in this case 
the analysis suggests that k > 1 + n does not occur unless dlQ^ vanishes 
asymptotically at the spike, for all integers < j < n. 

6 Further developments 

In section ^ we presented constructions of solutions of the Gowdy equations 
with spikes which are only local in x. Fortunately it is not hard to see that 
these spikes can be grafted onto solutions which are periodic in x and which 
therefore define global Gowdy spacetimes. Suppose that we have a local 
solution as constructed above with a true spike at x = x spike- Choose r 
large enough and \x — Xspikel small enough such that for a suitably chosen 
/c* the regions Ri and R2 defined in section ^ have the desired properties. 
Let 7i and 72 be the unique curves in the (r, x)-plane which are projections 
of null geodesies of the spacetime with y and z constant and which tend to 
X spike as T — > 00. Suppose r chosen large enough so that the parts of 71 
and 72 with values of r this large are contained in i?2- (To see that this 
is possible it is easier to think in terms of (t, x) coordinates, which make 
these curves into straight lines.) Choose some tq so that the straight line S 
given by r = To intersects Ri. Let 5*2 be the intersection of R2 with S. The 
restriction of the initial data on S for the given solution to a neighbourhood 
of 5*2 already determines the features of the solution which are characteristic 
of true spikes. The given data can be modified away from a neighbourhood 
of 5*2 and extended to data which are periodic in x. The extended data 
determine a global Gowdy spacetime which contains the original local spike. 
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In [|T^ Chrusciel obtained a sufficient criterion for the blow-up of the 
Kretschmann scalar near a point of the singularity in a Gowdy spacetime. 
This is his equation (3.4.5) which implies that the leading order behaviour 
of the curvature is as in the corresponding Kasner solution. It is interesting 
to see that this criterion is satisfied for the generic true spikes constructed 
above. This suggests that the techniques he uses to prove cosmic censorship 
in Gowdy spacetimes may admit extensions to wider classes, or even to the 
general case. Suppose that we have a solution of the Gowdy equations with 
a true spike at a; = of the kind constructed in section |^. Suppose further 
that the genericity condition Xqx{0) 7^ is satisfied. The criterion of |T2| , 
reexpressed in our variables, is 

/° /* (Pj + e^^Ql)dxdt < 00 (28) 

Jto J-t 

The region of integration is entirely contained in the region R2 when to is 
close to zero. Using the asymptotic expansion of the solution which is valid 
in that region shows that, up to terms of lower order, the integral to be 
computed is proportional to 

ft 

/ r^^x^dxdt (29) 

to J~t 

The integral is finite when k < 2, which is a condition which naturally arises 
from our construction. 

Next we produce spikes which have asymptotic velocity which is greater 
than that of the spikes produced in section ^, and which are consistent with 
the conditions suggested by the analysis in section |. An inductive argu- 
ment shows that successive applications of the transformations (0) and (0) 
produce spikes with \k\ arbitrarily large. To begin the inductive argument, 
let i be a nonnegative integer. Given a Gowdy solution labeled {P2i,Q2i), 

I GE 

let (P2i,Q2i) — ^ (^2i+i,Q2i+i) define (Ai+i, Q2i+i)- Let (P2i+i, Q2i+i) — ' 
{P2i+2;Q2i+2) define (F2i+2, (52i+2)- Then 

P2i+i = F2i + ln(Q| + e-2^-) 

= g2^ + ' -2P. (30) 



VTQ^.+l = -e''''''QlidrQ2i + drQ2i + 2idrP2i)Q 



2i- 
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And 



drQ2i+2 = e^^''''^-^~^Ql^,Q2^-e-^^{^,Q2^ + 2{^,P2,)Q2i) (31) 

dxQ2i+2 = G^'^^^'Qli 9TQ2i — drQ2i — 2 {drP2i)Q2i- 

Let (Poi Qo) be as in section ^ Then the following four properties hold for 
i = 0. 

1) lim dlQ2i = at Xgpike for all non-negative integers j, such that j < i. 

2) i < lim drP2i < i + ^ at Xs„ike, while < lim drP2i < 1 elsewhere in 
a neighbourhood of Xgpike- 

3) lim diiiQ^.^-^ = at Xspike for all nonnegative integers j such that j < i. 

4) -(z+l) < lim^^oo drP2i+i < -i at Xspike, while < lim^^oo drP2i+i < 1 
elsewhere in a neighbourhood of Xgpike- 

Now assume that these four properties hold for i equal to some nonneg- 
ative integer, n. Then from equations (|30|) and (0), and if the constant of 
integration for Q2n+2 is set such that limT-_,oo Q2n+2 = at Xgpike, it follows 
that the four properties hold for i = r?, + 1. So by induction, these four prop- 
erties hold for any nonnegative integer, i, with the constant of integration 
thus determined for each application of the GE transformation. 

The spikes just produced are consistent with the following conjecture, as is 
the analysis in section |^. Given some non-negative integer n, if —n > kgpike > 
— (n + 1), then at Xgpike, linir^oo f^^vrg = for all integers, < m < n. 
And given some positive integer n, if n < kspike < {n + 1) then at Xgpike, 
linv^oo 9^Q = for all integers, 1 < m < n. 

Note that {P2i+i,Q2i+i) is equivalent to {P2i,Q2i) since they are related 
by inversion. Thus, to determine whether the spikes just produced are true 
or false we need only examine the curvature of {P2i,Q2i)- If criterion (p8[ ) 
is satisfied, the curvature of (P2i, Q2i) near the singularity is as in the corre- 
sponding Kasner spacetimes. Therefore the discussion at the end of section ^ 
applies with < A; < 1 for x 7^ Xspike and i < k{xspike) < i + 1- It follows 
that (P„, Qn) represents a true spike if n > 1 and criterion ( pS] ) is satisfied. 

That criterion (pS]) is satisfied can be checked by obtaining expressions for 
{P2i,Q2i) on the region of integration. At Xgpike these are, for small enough 

P2i = {t + k)r + 0{l) Q2^ = X2. + o{e-'^) (32) 
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for some c > 0, and with di.X2i vanishing at Xgpike for all integers, 1 < j < i. 
These expressions follow from arguments similar to those concerning the 
expressions for {P2,Q2) on R2. The quantities which must be estimated on 
the region of integration are of the form e^"^~^'''''^\X2m\, with dlX2m vanishing 
at X spike for I < j < m. For x close enough but not equal to x spike, Pii = 
kr + 0{1) on the region of integration in criterion (|28|) for small enough to, 
with < A; < 1. The expression for Q2i is that given in (p2D. It follows that 
criterion (|28|) is satisfied in a neighbourhood of x spike- 

Because of the restrictive conditions which are satisfied at these higher 
velocity spikes, it is unlikely for them to have occurred in numerical simu- 
lations, since initial data for numerical simulations has not been specially 
chosen to obtain them. 

If we look beyond Gowdy spacetimes, the point of view of BKL and the 
method of consistent potentials suggest that there is a close association be- 
tween spikes and oscillations in time. According to the picture of BKL general 
solutions of the vacuum Einstein equations should show behaviour near the 
singularity similar to that of the mixmaster solution and, in particular, an 
infinite number of oscillations. This could lead to the generation of densely 
distributed inhomogeneous structure, which even calls the consistency of the 
BKL picture into question (cf. 0.) Issues of this type are beyond the reach 
of the rigorous mathematical techniques available at this time. However we 
can address one interesting issue. In it was proved, following a sugges- 
tion of Belinskii and Khalatnikov that coupling the Einstein equations 
to a scalar field can lead to mixmaster oscillations being suppressed. (See 
also the results on the homogeneous case in One might then ask if 

the presence of a scalar field could suppress the occurrence of (true) spikes. 
We will show that this is not the case. In fact this is apparent within the 
BKL approach. The scalar field forbids the occurrence of an infinite number 
of oscillations and hence the occurrence of an unlimited number of spikes. 
However it does not forbid finitely many oscillations and finitely many spikes. 

The occurrence of true spikes in solutions of the Einstein-scalar field sys- 
tem will now be shown rigorously. The idea is to take a solution of the 
Gowdy equations with a true spike and produce from it a solution of the 
Einstein-scalar field equations with Gowdy symmetry which also contains a 
spike. The possibility of doing this is based on the structure of the field 
equations in this case. The metric can be written in the same form as in the 
vacuum case. The field equations look very similar to the equations in the 
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vacuum case. The equations for P and Q are unchanged. This is because 
they come from the projection of the Ricci tensor onto the orbits of the group 
action. Since the Ricci tensor is proportional to the tensor product of the 
gradient of the scalar field with itself, and this gradient has vanishing projec- 
tion onto the orbits, the corresponding components of the Einstein equations 
are unchanged. The scalar field satisfies the equation 

4>TT = e~'^''(j)a:x (33) 

Note that this is the same equation as satisfied by P in the case that Q 
is zero. Thus there is a strong connection to polarized Gowdy solutions. 
In particular, we understand the asymptotics of solutions of the equation 
for (f) completely. The equations for A are modified by adding to the right 
hand sides contributions from (p formally identical to those due to P. Let 
us start with a solution of the Gowdy equation with a true spike and take 
any non-vanishing solution of the equation for 0. We can then construct a 
solution of the Einstein-scalar field equations as just indicated. The function 
P has not been changed and so still has a spike. The calculation for the 
curvature is essentially the same as in the vacuum case. In the calculation 
for the general case the leading order is given by the Bianchi I case. The 
function is asymptotically of the form Alnt + B. The only difference in 
the computation of the Kretschmann scalar is that the Kasner exponents 
are replaced by the generalized Kasner exponents satisfying the equation 
Pi + pi + pI — 1 ~ IGnA^. If A is small then the limits of the generalized 
Kasner exponents of the solutions with scalar field at the singularity will differ 
only a little from the limits of the Kasner exponents of the original vacuum 
solution. It follows that the addition of a scalar field does not always remove 
the discontinuous limiting behaviour of the Kretschmann scalar. 

The fact that we were able to obtain so much information about spikes 
in this paper was based on a trick (the use of the inversion and Gowdy-to- 
Ernst transformation). Spikes can be expected to occur in much more general 
spacetimes but this trick cannot be expected to extend very far. The method 
of consistent potentials is not limited in the same way and it would be very 
desirable to transform that technique into a tool which can be used to prove 
theorems. Having one case, namely the Gowdy case, which is understood in 
detail may provide valuable guidance as to how to do this. 
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